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In this work, we apply black box methods (methods not requiring input) to find 
excited-state energies. A variety of such methods for lattice QCD were introduced 
at the 3rd iteration of the numerical workshop series. We first review a selection of 
approaches that have been used in lattice calculations to determine multiple energy 
states: multiple correlator fits, the variational method and Bayesian fitting. In the 
second half, we will focus on a black box method, the multi-effective mass. We 
demonstrate the approach on anisotropic lattice data, extracting multiple states 
from single correlators. Without complicated operator construction or specialized 
fitting programs, the black box method shows good consistency with the traditional 
approaches. 
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I. INTRODUCTION 

Lattice QCD has successfully provided many experimental quantities from first-principles 
calculation, sometimes predicting experimental results before they are measured. However, 
its success has mostly been restricted to measurements of ground-state quantities. In spec- 
troscopy, for example, there are many poorly known states which require theoretical input 
to be identified. At the EBAC at Jefferson Lab, dynamical reaction models have been 
developed to interpret extracted N* parameters in terms of QCD HQ. 

Among these excited nucleon states, the nature of the Roper resonance, iV(1440) or N', 
has been the subject of interest since its discovery in the 1960s. It is quite surprising that 
the nucleon's excited-state mass is lower than its negative-parity partner, a phenomenon 
never observed in meson systems. There are several beyond-the-quark model interpretations 
of the Roper state, for example, as the hybrid state that couples predominantly to QCD 
currents with some gluonic contribution 3J or as a five- quar k (meson-baryon) state [4]. Earlier 
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ll| . found a spectrum inverted with 



respect to experiment, with the P\\ N' heavier than the opposite-parity state Su. However, 
this has been clarified by the study in Ref. jf| with a larger lattice box « (3.2 fm) 3 and 
a broad pion mass study from 870 MeV to 180 MeV. They found that the mass of the 
N' becomes lighter than the Su when the pion mass gets lighter than 300-400 MeV. The 
findings of Ref. [(J indicate a rapid crossover of the first positive- and negative-parity excited 
nucleon states close to the chiral limit. We need a full dynamical lattice QCD calculation 
with proper extrapolation to (or calculation nearby) the physical pion mass to resolve the 
Roper issue once and for all. 

The excited transition form factor is also an interesting quantity for lattice QCD to 
explore. Experiments at Jefferson Laboratory, MIT-Bates, LEGS, Mainz, Bonn, GRAAL, 
and Spring-8 offer new opportunities to understand in detail how nucleon resonance (N*) 
properties emerge from non-perturbative aspects of QCD. Consider the Roper transition 
form factor, for example. Figure [I] shows preliminary data from CLAS collaboration that 
covers a large range of photon virtuality Q 2 . In the region Q 2 < 1.5 GeV 2 , both the 
transverse amplitude Ai/ 2 (Q 2 ), and the longitudinal amplitude Si/ 2 (Q 2 ), drop rapidly in 
magnitude. This is well described in relativistic quark models with light-cone dynamics, 
and the sign is consistent with the non-relativistic version. However, in the low-Q 2 region, 
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FIG. 1: Preliminary results on the transverse (left) and longitudinal (right) amplitudes of the 
Roper resonance. The solid lines represent non-relativistic quark model predictions, the dashed 
lines are predictions of relativistic quark models. 



Ai/ 2 (Q 2 ) becomes negative; this is not understood within constituent quark models and 
requires inclusion of meson degrees of freedom. There are various QCD-based hadron models 
such as the well developed constituent quark model[12|, |l3(] and the covariant model based on 
Dyson- Schwinger Equations |l4j. A model-independent study of these quantities from lattice 
QCD will serve as valuable help to phenomenologists in analyzing experimental data from 
stronger theoretical ground to gain a better understanding of low-Q 2 physics. Recently, we 
have submitted a proposal for studying this physics to the USQCD collaboration; we expect 
the first lattice QCD calculation of this quantity fairly soon. 

Even when one only cares about getting ground-state quantities, one cannot quote precise 
numbers without estimating excited-state contamination. By explicitly adding excited-state 
degrees of freedom, we can obtain more accurate assessment of ground-state quantities. 
Examples will be given in later sections. 

The structure of this proceeding is as follows: Sec. [TTlgives a brief lattice QCD introduction 
for the non-latticists. We discuss various approaches that have been used in lattice QCD 



for extracting excited states, such as mu 
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tiple-correlator least-squares fitting, the variational 
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191 ]. and last but not least, black box methods[20], 
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including both analytic excited effective masses and linear predictions. Conclusions and 
future outlook are given in Sec. [IVl 

II. LATTICE QCD 

This section provides a brief introduction to lattice QCD calculations for non-experts. 

Quantum chromodynamics (QCD), a non-Abelian (Yang-Mills) gauge theory with gauge 
group SU(3), is currently the best description of the strong interaction. In this theory, the 
physical observables are calculated from the path integral 

(n\o\n) = \ l[du][d#^e- s ^ u ^- SG ^o{u^^ (i) 

Zj J 

in Euclidean spacetime. The asymptotic freedom of QCD is a consequence of the negativity 
of the beta function for small numbers of flavors or large numbers of colors in non-Abelian 
gauge theory. When the coupling becomes weak at short distances, we can perform pertur- 
bative calculations to extract physics analytically. However, at large distances the strong 
coupling which accounts for quark confinement leads to a perturbative series that no longer 
converges, and we need a non-perturbative treatment, like lattice QCD, to extract physics 
from QCD. 

Lattice QCD is a discrete version of continuum QCD theory. The path integral over 
field strengths at infinitely many spacetime points from continuum QCD is approximated 
using only finitely many points in a discrete Euclidian space lattice with periodic boundary 
conditions. In momentum-space, this automatically provides both an ultraviolet cutoff at 
the inverse lattice spacing and an infrared cutoff at the inverse box size. 

By discretizing spacetime, we can then use Monte Carlo integration combined with the 
"importance sampling" technique to calculate the path integral numerically. At the end of 
the day, we need to take a — > and V — > oo to reach the continuum limit. However, to 
simulate at the real pion mass (while at the same time keeping the lattice box big enough 
to avoid finite- volume effects) would require much faster supercomputers than are currently 
available to the community. Thus, we normally calculate with a few unrealistically large 
values of the pion mass and then use low-energy chiral effective theory to extrapolate to the 
physical pion mass. 

Most observables calculated in lattice QCD can be categorized as two-point and three- 
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point Green functions: 

C 2 r pt (X src ,X) = J2^(J(X)J(X src )) a ^ (2) 

a,/3 

Cfptl^src, X, X sn k) = T a,/3 ( J (X sn k) O (X) J (X src ) ) a> p (3) 

After taking momentum (and spin for baryons) projection (ignoring the variety of boundary 
condition choices possible), these become 

C 2pt (t) = Y,Up)e- E "\ (4) 

n 

C^t U t,t f ,%p f ) = X)E^'(P/)^(Pi) e " ( * / "* ,B;,( ^ )e " (t ~ ti) ^ (?<) X fOTm faCtOTS > ( 5 ) 

n n ' 

where Z contains not only the overlap factor between the vacuum and the operator, but 
also the kinetic factor with energies or masses. Note that in either of the two- or three-point 
functions, the different states are mixed, and the signal decays exponentially as a function 
of t. This gives us a challenge: we want to disentangle the different channels, and this must 
be done before the signal dies out at large times. 

One usually obtains ground-state observables at large t, after the excited states die out, 
a method which requires large time dimension. The lighter the particle is, the longer the t 
required to get a "good enough" ground state. As lattice simulations go to light sea pion 
masses and finer lattice spacings, even getting an accurate ground state would not be easy 
without including excited states in the analysis. 

Moreover, sometimes the large-t solution does not work well, since the signal-to-noise ratio 
rapidly decreases. For example, when LHPC & SESAm|2i| calculated the quark helicity 
distribution, they found a 50% increase in error at t scp = 14 over the error at separation 12. 
Confronting the excited states might be a better strategy than avoiding them, especially for 
precision calculations in the future. 



III. PROBING EXCITED STATES 



There are various approaches that have been used widely by lattice theorists to probe 
excited states: multiple-correlator fits, variational methods and Bayesian fitting methods. 
We would particularly like to continue the discussion of black box methods and show how 
they may be implemented on lattice QCD data for the first time. 
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FIG. 2: Effective mass plots from the 7 smeared-point proton correlators used in this work 

Note that the data on which we demonstrate these methods is from a "quenched" (that 
is, ignoring fermion loops in the sea) study on 16 3 x 64 anisotropic lattices with Wilson 
gauge action and nonperturbative clover fermion action. The spatial lattice spacing a s is 
aboutO.125 fm with anisotropy 3 (that is, temporal spacing a^ 1 « 6 GeV). Specifically, our 
data are proton correlators using 7 gaussian smearing parameters, 0.5-6.5 in steps of 1.0, 
including both smeared-point and smeared-smeared source-sink combinations. 



A. Multiple-Correlator Fits 

The brute-force way of extracting excited states is to fit using a form that includes them 
explicitly. Usually, one needs multiple correlators as input to make the fit program converge. 
With multiple smeared correlators, one minimizes the quantity 

X = ? m ' (6) 

where the correlators G s are each fit by the sum of n exponentials, having independent 
amplitudes but the same masses (energies) E. To extract n states, one at needs at least n 
distinguishable input correlators. 

Using this technique, we find E = 0.240(3), E l = 0.40(3) and E 2 = 0.67(4). 



B. Variational Method 



We construct an r x r spectrum correlation matrix, Cij(t), where each element of the 
matrix is a correlator composed from different smeared sources or operators Oi and Oj. 
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Then we consider the generalized eigenvalue problem 

C{t)iP = \(t,t )C(t )il>, (7) 

where to determines the range of validity of our extraction of the lowest r eigenstates. If 
to is too large, the highest-lying states will have exponentially decreased too far to have 
good signal-to-noise ratio; if to is too small, many states above the r we can determine will 
contaminate our extraction. Over some intermediate range in to, we should find consistent 
results. 

If the eigenvector for this system is |a), and a goes from 1 to r. Thus the correlation 
matrix can be approximated as 

= £ tT«i e-°S (8) 

n=l 

with modified eigenvalue 

\ n (t,t ) = e- (t - to)En , (9) 

by solving 

C(t )- 1/2 C(t)C(t )- 1/ ^ = A(t, t )V- (10) 

Here we apply the above steps to a 7 x 7 smeared-smeared proton correlator matrix. 
We solve the eigensystem for individual A n , and use these eigenvalues to construct effective 
mass plots, as shown in Figure [31 The masses of the lowest four eigenstates as fit in the 
regions shown in the figure are: E = 0.2526(24), E l = 0.406(7), E 2 = 0.533(16) and 
E 3 = 0.633(16). 

C. Bayesian Fitting 

One method for extracting additional masses from a single correlator is to apply more 
sophisticated fitting techniques. Using information we already know, we can introduce what 
are known as Bayesian priors to constrain our fits. These constraints may improve the 
precision with which difficult-to-determine masses are found by a fitting method. 

The Bayesian technique applies known information to shift the estimated probabilities of 
correctness for parameter sets being considered by a fitting procedure. This prior information 
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FIG. 3: Effective mass plots from generalized eigenvalues of a 7 x 7 smeared-smeared proton 
correlator matrix from anisotropic lattices. 

may be as simple as the assertion that all physical masses must be positive (Ei > 0) or that 
each new excited state must be heavier than the states preceding it (Ei > In our 

case, we will attempt to extract states from a correlator in an iterative fashion. Assuming 
we have somehow acquired a ground state for our data (say, by a non-Bayesian fit), we can 
make the prior assumption that the ground state is actually equal to that fitted value within 
the determined uncertainty. 

Using this prior information, we wish to modify the likelihood estimate for all parame- 
ter sets considered by our Bayesian fitter. The probability of correctness for a particular 
parameter set is generally represented by the x 2 value, which for a non-Bayesian fit is 



Xdof 



1 



N -n 



N 

E 



Vi - f({a,E}o..n-i,Xi) 



(11) 



i=x \ 6 Vi J 

where we fit N data points y ± 8y to some form / parametrized by n amplitudes A and 
masses (energies) E. Rather than allowing the fitter to arbitrarily change the amplitudes and 
masses, we would like it to take into consideration the masses we have already determined. 
Therefore, we simply add a term that will increase x 2 (decrease the probability of correctness) 
if the fitter moves away from these values: 

A 2 
m— 1 / rp t? \ 

a * 2 =5(t)- (12) 

where we have already determined the first m masses to be E ± 5E. 

We apply such a Bayesian fit to the most pointlike of the correlators mentioned in the 
previous section. The results are shown in Figure HI In the top-most plot, an ordinary fit 
determines the ground-state mass E = 0.243(5) over a short plateau t G [35,47]. Using 
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FIG. 4: A series of effective mass plots showing (red) the standard ground-state fit, (blue) an 
excited-state fit using the ground state as a Bayesian prior, and (green) the second-excited fit 
using ground and first-excited masses as priors. 

this prior information, we extend the fit region to t € [25,47], finding Eq = 0.2421(22) and 
Ei = 0.47(11). Finally, we apply both these masses to the fit region t G [20,47], finding 
E = 0.2420(8), Ex = 0.463(26) and E 2 = 0.7(9). Since our second-excited mass is not 
statistically distinguishable from noise, it does not make sense to continue this procedure to 
yet-higher masses. 
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D. Black Box Method 



In the 3rd iteration of this workshop, G. T. Fleming[20(] talked about "black box methods" 
used in NMR (nuclear magnetic resonance) imaging, where practitioners deal with problems 
of the form 

K 



+ e n . 



k=l 



This is similar to lattice correlators, which have the general form 

K K 
n=l n=l 

where a n = e~ En . This forms a Vandermonde system 

/ </>i(ti,a) • ■ • (/> K (ti,a) ^ 



y = $a, ; where $ 



\(f)x(t N ,a) ■ ■ ■ <f) K (t N ,a) J 



(13) 



(14) 



(15) 
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Given a single correlator (with sufficient length of t) , one can in principle solve for multiple 

a's and a's 



2/2 



\ 



N N 



a 2 



(16) 



The question is: can one solve for N/2 states from one correlator of length N? 



1. Analytic approach 
In the K = 1 case, the system dominated by the ground state, 

a i 



( V, ^ 




u ) 


K Vt+i ) 


- 





(17) 



The solution is easy a% = yt+i/yt- This is familiar to lattice theorists as the "effective mass" 
M e fj = \og(y t +i/yt)- Applying this formula at various times produces an effective mass 
plot, which is helpful in determining appropriate fitting ranges, away from excited-state 
contamination. 

Two states (the K = 2 case) may be extracted using 



fj2 
V3 



4 4 
\ «1 «2 j 



X 




(18) 



which leads to a more complicated solution 



"1,2 



2/12/4 - 2/22/3 ± V (2/22/3 - 2/12/4) 2 + 4(y| - 2/12/3) (2/22/4 



2/1) 



2(2/i2/3 



2/1) 



(19) 



One thing to note that is that the modified correlator approach, proposed by D. Gaudag- 
nol et al. [Tj] is mathematically identical to the product of our solutions 0^2, giving the sum 
E + E,. 

An example using our nucleon data is given in Fig. 

Such extensions of the effective mass technique provide simple but slightly different points 
of view from what we normally do in lattice calculations. Simple algebra gives us multiple 
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FIG. 5: Higher-effective mass plots with (top-to-bottom) one, two, three and four masses. The 
colors indicate the (black) ground, (red) first-excited, (blue) second-excited and (green) third- 
excited states. 

states from a single correlator, which might be difficult to achieve using least-squares fitting. 
It is possible to calculate the second-excited state in similar fashion. However, there is 
a limitation to this analytic approach: Abel's Impossibility Theorem, which states that 
algebraic solutions are only possible for N < 5. To obtain more states, we would need to 
solve the resulting polynomial equations numerically. 
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2. Linear prediction 



If we are interested in extracting K states, we can first construct a polynomial with 
coefficients 



n n=1 (a-a 



K 

n > ~= Pi 
n=l 



K-n 



a 



(20) 



We can "predict" a data point at a later if in terms of the earlier time points: 



K 



y t i — — PnOt 1 n where m > t' 



(21) 



n=l 



for ideal data. This is called a "linear prediction". However, for non-toy data, background 
fluctuations will be present; thus, Eq. [2T1 becomes 

M 



YsPnO?' n + V t > 



(22) 



n=l 



Therefore, to extract the multiple states, we need to solve the Vandermonde system 



VM+l 



( 



Vo ■■■ Vm-i 

yi ■■■ yn 



\ ( „ \ 

Pm 
Pm-i 



where N > 2M. 



(23) 



\ VN-l J \ Vn-m-i ■ ■ ■ VN-2 J \ Px J 

We demonstrate this approach using the smallest smearing parameter (0.5) smeared-point 
correlator. There are a few parameters in the linear prediction method which we can tune: 
the number of desired states K, the number of time slices used to predict the later time point 
N, and the order of the polynomial M. In this work, we will show a selection of the better 
choices in these degrees of freedom. Figure M shows the effective mass plot for K = 2,3,4 
from a single Gaussian smeared-point corrleator with fixed parameters iV = 20 and M — 8. 
The excited states are consistent with each other as one increases the value of K. Since we 
have used a large value of N to form the polynomial, each point uses information extracted 
from 20 time slices. Thus, one does not need a large plateau to determine the final mass. 
One also notes that since we only use a single correlator to extract multiple states, the 
multiple states will be correlated; that is, large errors on higher-excited states will make 
the ground state noisy as well. A future improvement would naturally be to extend this 
approach to multiple correlators. 
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We compare this linear prediction approach with the variational method (4 x 4, with 
smearing parameter ranging 0.5-3.5), as shown in Figured Here we shift time with respect 
to the linear prediction plot by 10 to have a better comparison with the plateau region from 
the variational method. For the ground state, the numbers are consistent with the result 
from the varitional approach, including the size of the errorbar. This is remarkable, given 
that the amount of input information is a factor of 16 less in the linear prediction approach. 
The first-excited state is consistent but has larger errorbar, which is no surprise. As for the 
second-excited state, it seems to be consistent with the varitional ones but definitely needs 
more stastictics. The third-excited state is much larger than expected from the varitional 
approach, which might be caused by contamination from even higher excited states. 

IV. CONCLUSION AND OUTLOOK 

There is a lot of interesting physics involving excited states, but they are traditionally 
difficult to handle in lattice QCD. Even for physicists who are only interested in the ground 
states, the estimation of excited states can improve the precision of their calculations. 

In this work, we adopt NMR-inspired methods which provide an interesting alternate 
point of view for looking at lattice QCD spectroscopy. We successfully extract multiple 
states from single correlators, working with quenched anisotropic lattice data. We also 
demonstrate various other methods for determining excited states (multiple-correlator fit, 
variational method and Bayesian fit). These methods require different amounts of input 
information and manual guidance, which may introduce systematics that are difficult to 
quantify. The black box methods provide alternatives without any ambigiuity and without 
the need for careful tuning. 

Future work will extend this approach to multiple correlators [22], introducing more in- 
formation to improve the signal. 
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FIG. 6: Effective mass plots from the linear prediction black-box method at fixed parameters 
N = 20 and M = 8 



performed on clusters at Jefferson Laboratory under the SciDAC initiative, using the Chroma 
software suite [231] on clusters at Jefferson Laboratory using time awarded under the SciDAC 
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FIG. 7: Comparison between the variational method (4 x 4) and linear prediction with 4 extracted 
states 
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